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Abstract 

The differences between lattice Monte Carlo data and perturbation 
theory are usually associated with the "bad" behaviour of the bare lattice 
coupling g$ due to the effects of large (and unknown) higher order coeffi- 
cients in the go perturbative series. In this philosophy a new, renormalised 
coupling g' is defined with the aim of reducing the higher order coefficients 
of the perturbative series in g' . An improvement in the agreement between 
lattice data and this new perturbation series is generally observed. 

In this paper an alternative scenario is discussed where lattice artifacts 
are proposed as the cause of the disagreement between lattice data and the 
go-perturbative series. We find that this interpretation provides excellent 
agreement between lattice data and perturbation theory in go corrected 
for lattice artifacts. We show that this viewpoint leads typically to an 
order of magnitude improvement in the agreement between lattice data 
and perturbation theory, compared to typical g' perturbation expansions. 
The success of this procedure leads to a determination of ° of 220 ±20 
MeV. Lattice data studied includes quenched values of the string tension, 
the hadronic scale ro, the discrete beta function A/3(/3), M p , f n and the 
IP — IS splitting in charmonium. The new 3-loop term of the lattice 
/3-function [Q] has been incorporated in this study. 

A discussion of the implication of this result for lattice calculations is 
presented. 
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1 Introduction 



A necessary condition for lattice predictions of QCD and other asymptotically 
free theories to have physical (continuum) relevance is that they reproduce weak 
coupling perturbation theory (PT) in the limit of the bare coupling go — > 0. 
This perturbative scaling behaviour (a.k.a. asymptotic scaling) has not yet been 
observed for complicated theories like QCD. Even for simple systems such as 
the 2-dimensional 0(3) model, there is evidence that in present simulations, the 
correlation length, £, is still 4% away from the perturbative behaviour even for 
huge values of f ~ 0(H) 5 ) §. 

These comments apply when the "naive" , bare lattice coupling from the lat- 
tice action, g , is used as the expansion parameter in the perturbative series. 
As a result of this disappointing disagreement, various workers have proposed 
methods of improving the convergence of the perturbation series by re-expanding 
it in terms of some new coupling g' []3], || . They argue that the lack of (2-loop) 
perturbative scaling is due to the presence of higher order terms in the pertur- 
bation expansion appearing with large coefficients. If a more physical coupling, 
g', could be found, then, it is argued that the perturbation series expressed in 
g' will converge faster (i.e. it will have smaller higher order terms). There has 
been evidence to support this philosophy (see for example 0, §J)- We refer to 
this method as the "renormalised coupling" approach. 

This paper studies an alternative viewpoint in which the above disagreement 
stems from lattice artifacts, i.e. cut-off effects 0. All lattice Monte Carlo data is 
obtained on a lattice with finite lattice spacing, a, with an action that is correct 
to 0(a n ) for some n. Therefore, a priori, corrections of the form of 0(a n ) at best, 
are present in any observable. These terms are usually treated in the final stage 
of the analysis, when an extrapolation to the continuum limit of some physical 
quantity is performed (see for example |5|, [7|, ||). In this paper, these terms are 
shown to provide the mismatch between lattice Monte Carlo data, expressed in 
lattice units, and PT without resorting to the use of a renormalised coupling g' . 
That is, bare lattice data is reproduced by go-PT with the simple addition of 
terms of 0(a n ), with an appropriate value for n. We call this agreement between 
lattice data and PT "Lattice-Distorted Perturbative Scaling" . 

The QCD quantities studied in this analysis are: the string tension, y/a; the 
hadronic scale, r |§; M p ; f n ; the IP — 15 splitting in charmonium |7|; and 
the discrete beta function A/3. We find that the lattice-distorted PT reproduces 
the Monte Carlo data for all the quantities considered. When fits are performed 
to x/a, r and A/3 using various proposed renormalised coupling schemes, the 
X 2 /dof are around an order of magnitude worse than those obtained with lattice- 
distorted PT. [] The implication of this result is that the higher order terms in 

1 In the case of M p , / w and the IP — 15 splitting, the lattice data is too noisy to constrain 
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the <?o-PT expansion of these quantities do not "contaminate" the data, whereas 
the (D(a)-type terms do. 

As this paper was in the final stages of preparation, a calculation of the 3-loop 
coefficient of the lattice /5-function appeared Jl|. We has included this term in 
this study where appropriate. 

The plan of this paper is as follows. In the next section, the (straightforward) 
formalism of lattice-distorted PT is presented. Section || fits the QCD lattice data 
to lattice-distorted PT and obtains a value of = 220 ± 20 MeV. Section 

|4] fits the same data to the renormalised coupling-PT schemes, specifically the 
9ms scheme of [J7), the gv-PT scheme of [§J, and the g^-PT schemes of & 0. 
We conclude with an interpretation of these results, and comment on how this 
method can be used to find the continuum value of quantities determined on the 
lattice. 

A brief account of the ideas in this paper appears in [[TTJ. 



2 Lattice-Distorted Perturbation Theory 

We begin this discussion with one of the fundamental quantities in perturbative 
field theory, the (3-loop) beta function: 

do 2 

P(9 2 ) = -a=jfc = -2b g 4 - 2b l9 & - 2b 2 g\ 

where the one- and two-loop coefficients for the quenched theory are 

11 , 102 



(4tt) 2 ' 1 (4vr) 4 ' 

The (scheme- dependent) three-loop coefficient has recently been calculated [|TJ 
for the SU(N) lattice action, 

( N W Q ««o , 1433 - 8 2143 - 



2 V16tt 2 / \ ' N 2 N 4 

The /5-function can be intergrated to give the usual form for the running of the 
coupling with the U.V. cut-off 1/ 



a. 



a' 1 = (1) 



wherej 2 ] 



the fits. 

-h 

2 Note that the definition of fpx differs from that in |(| and |T^] by a factor of b to 
conform with convention. 
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f PT (g*) = e 55* (6 /)4 (l + d 25 2 ), (2) 

where c?2 = oZ3-(&i — &2&o)- In the case of the lattice scheme, we have d\ = 
J^(b{-b%) = 0.1896 1. 

All coupling and fields on the lattice, and therefore all observables, are di- 
mensionless. Lattice calculations set the scale a -1 by calculating some quantity 
on the lattice, eg. the string tension, a^, and comparing it with its experimental 
(dimensionful) value, <J exp : 



(J, 



exp 



<<7L 

Typical values for a~ l from recent simulations are listed in the third column of 
table |l] together the references and corresponding j3 = Q/g^ values in columns 
1 and 2. It is now easy to check if a" 1 follows go-PT (i.e. eq.(fj) with g — go)- 
Plotting A CT = oi^ 1 f pt{9q) * n fig IE we observe a non-constant behaviour, signaling 
the failure of Monte Carlo data to follow even 3-loop go-PT. (The failure of 2-loop 
perturbative running has already long been noted.) 

There are a number of possible causes of the disagreement. 



• quenching 

• finite volume effects 

• unphysically large value of the quark mass (relevant for "hadronic" quan- 
tities such as masses and decay constants, rather than for a as used in the 
above example) 

• a real non-perturbative effect 

• the inclusion of only a finite number of terms in the PT expansion 

• lattice artifacts due to the finiteness of a 



For the reasons outlined in ||, the first three effects cannot give rise to 
the sizeable discrepancy between lattice data and PT in fig |I|. (For example, 
quenching should modify A by an overall constant factor.) As far as true (i.e. 
continuum) non-perturbative effects are concerned, the overwhelming expectation 
is that for cut-offs of a -1 ^ 2 GeV these effects should be minimal. (In any case 
these effects have the same form as lattice artifacts since they are of the form 
e -1 /^ ~ 0{a).) Therefore, we can assume that the disagreement is due to either 
or both of the last two possibilities. 
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The effects of higher orders in go and the finiteness of a can be parametrised 
in the lattice beta function as follows: 



PM) = -a L j^ = -(2b g* + 2^ 6 + 2« + E b ti9o +2 ) 

x (l + E c nG?oK(<7o)) (3) 

n=l 

Here the 63 , 64 , ... are the (unknown) higher loop coefficients of the lattice beta 
function. (They are presumed to be large in the renormalised coupling approach.) 
The c n are the (non-universal) coefficients of the 0(a n ) pieces and are, in general, 
polynomial functions of <7q.[] 

Eq-dD can be integrated giving 

^(So) = 7%T x (1 + E dlA^Y 1 x (1 + E <M)fM)) (4) 

JPT{9o) i =A n=l 

Note that the 0(a n ) term has been expressed in terms of fp T . This can be done 
without any loss of generality since any difference between and fpx is higher 
order and can be absorbed into the coefficients c' m (gQ) for m > n. 

In the following section we study lattice-distorted PT by setting the higher 
order coefficients, d L , to zero, and fitting the data in order to determine the c' n . 
We perform this fit for both the 2-loop function (i.e. we set d\ to zero in /pt) 
and the 3-loop function. In section | the renormalised coupling ideas are studied: 
i.e. all the c' n are set to zero, and g is replaced by gj^ 0, gy 0, gE an d gE2 
0. In these cases, we also fit to the appropriate 2-loop and 3-loop formulae. 



3 Fits Using Lattice-Distorted PT 

In this section we fit lattice Monte Carlo data for a, r ||, M p , f v and the IP— IS* 
splitting to the form eq.(^) with the coefficients d L set to zero appropriately (i.e. 
we ignore higher order effects in g^). The data for these quantities is displayed 
in table [I] together with the references in column 1. 

We perform two fits to eq.(^). The first with fpj- defined with d\_ x = for 
/ > 3 (i.e. including only 2-loop terms - see eq.(||)). The second is with f PT 
defined with d\_ x = for / > 4 (i.e. including the newly calculated 3-loop term, 
b L 2 i). 

With the exception of the IP — IS splitting J7| which comes from the clover 



action [12, 13], we consider results only from the Wilson action [11]. This is 



We have used the replacement g% log a ~ 1 in eq. (||) since the difference between g$ and 



log 1 a can be incorporated into the higher order terms. 
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"Rel f3 = Q/g 2 a a l aJ aj aT 1 ~ T 



*Tg M p "/j- "IP-IS 



T5[ - WUP 5.5 0.80(1) 



T§ - WUP 5.6 0.91(2) 



T§ - MT C 5.7 1.073(3) 



1§ - WUP 5.7 1.14(2) 



Tjj - GF11 5.7 1.42(2) 1.08(5) 



7j - FNAL 5.7 1.15(8) 
m - WPC 5.74 1.44(3) 



T| - MT C 5.8 1.333(6) 



T§ - WUP 5.8 1.45(2) 



16 1 - MT C 5.9 1.63(2) 



T§ - WUP 5.9 1.85(5) 



7| - FNAL 5.9 1.78(9) 

17| - GF11 5.93 1.99(4) 1.78(5) 

g - WPC 6.0 2.25(10) 

| - WUP 6.0 1.94(5) 

g - UKQCD 6.0 2.04(2) 



T5| - WUP 6.0 2.11(2) 

g - UKQCD 6.0 2.19(4) 

2g - APE 6.0 2.23(5) 1.98(8) 

211 " APE 6 -° 2.18(9) 1.76(8) 

22[ - APE 6.1 2.64(16) 2.57(15) 

7j - FNAL 6.1 2.43(15) 

T7| - GF11 6.17 2.77(4) 2.56(7) 



| - WUP 6.2 2.72(3) 
g - UKQCD 6.2 2.73(3) 



T5| - WUP 6.2 2.94(2) 

1| - UKQCD 6.2 2.92(6) 

5g - APE 6.2 2.88(24) 2.69(24) 

- WPC 6.26 3.69(32) 

- WUP 6.4 3.62(4) 



T§ - WUP 6.4 3.95(3) 



53 - ELC 6.4 3.70(15) 3.7(6) 



22 1 - APE 6.4 4.09(18) 3.48(16) 



g - UKQCD 6.4 3.62(7) 



1| - UKQCD 6.4 3.90(6) 



2J - UKQCD 6.5 4.12(4) 



| - WUP 6.8 6.0(1) 

§ - WUP 6.8 6.7(2) 



Table 1: Values for a 1 obtained from various group's work using the Wilson 
action. Note that we have used a exp = (440 MeV) 2 , and r^ 1 = 400 MeV. 
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a 1 from 


Af3((3) 






M p 


u 


IP - 15 


<7 -Perturbation Theory 


A L [MeV] 
X 2 /dof 


3.856(8) 
484 


4.91(2) 
262 


5.00(4) 
10 


4.54(6) 
9 


4.5(2) 
6 


702 


Leading-Order Lattice Distorted PT 

i.e. Fit using eqs.(|||7|) 


A L [MeV] 
y 

X 2 /dof 


5.85(3) 

0.204(2) 

3 


6.02(3) 

0.150(2) 

16 


6.6(2) 

0.22(2) 

1.1 


6.9(3) 

0.34(3) 

1.6 


7.6(9) 

0.35(6) 

0.3 


0.195(2) 
3.5 


Next-to-Leading-Order Lattice Distorted PT 

i.e. Fit using eqs.(|].[7]) 


A L [MeV] 
Y 

X 2 /dof 


6.02(5) 
0.26(2) 
-0.024(6) 
1.7 


6.58(5) 
0.29(1) 
-0.046(3) 
1.4 








0.252(6) 

-0.025(3) 

1.7 


{/^^-Perturbation Theory 


Am^ [MeV] 
X 2 /dof 


53.3(1) 
160 


64.2(2) 
47 


65.7(5) 
1.3 


64.5(8) 
2.5 


59(2) 
1.5 


78 


^^-Perturbation Theory 


M v [MeV] 
X 2 /dof 


81.2(2) 
176 


98.2(3) 
54 


100.3(7) 
1.4 


99(1) 
2.7 


91(3) 
1.6 


89 


<7^^ -Perturbation Theory 


M v [MeV] 
X 2 /dof 


104.3(2) 
31 


118.9(4) 
13 


123.9(9) 
5.2 


119(2) 
1.4 


113(4) 
0.13 
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<7£-Perturbation Theory 


A E [MeV] 
X 2 /dof 


14.80(3) 
52 


17.09(6) 
15 


17.7(1) 
3.6 


17.1(2) 
1.4 


16.1(6) 
0.3 


19 


<7.E2-Perturbation Theory 


A E [MeV] 
X 2 /dof 


8.07(2) 
39 


9.25(3) 
12 


9.62(7) 
4.4 


9.3(1) 
1.4 


8.8(3) 
0.2 


13 



Table 2: Fits of lattice data to 2-loop PT using i) g^-PT, ii) Leading Order 
Lattice Distorted PT, Hi) Next-to-Leading Order Lattice Distorted PT, iv) gjj s - 
PT, v) gy' H ^-PT, vi) gE-PT and vii) gsp-PT. 





a 1 from 


A/3(/3) 








u 


IP - IS 


<7 -Perturbation Theory 


A L [MeV] 
X 2 /dof 


4.62(1) 
448 


5.85(2) 
239 


5.96(4) 
9 


5.40(8) 
8 


5.3(2) 
5 


625 


Leading-Order Lattice Distorted PT 

i.e. Fit using eqs.(||[3) 


A L [MeV] 
X 2 /dof 


6.86(4) 

0.193(2) 

3 


7.07(3) 

0.141(2) 

15 


7.7(2) 

0.20(2) 

1.1 


8.0(4) 

0.32(3) 

1.6 


8(1) 

0.34(6) 

0.3 


0.184(2) 
3 


Next-to-Leading-Order Lat 

i.e. Fit using eqs 


tice Distorted PT 
•d® 


Al [MeV] 
Y 

TIM 

X 2 /dof 


7.01(6) 
0.24(1) 
-0.019(6) 
1.8 


7.68(6) 
0.27(1) 
-0.042(3) 
1.4 








0.23(1) 

-0.021(4) 

1.6 


gf^-Perturbation Theory 


A E [MeV] 
X 2 /dof 


17.02(4) 
36 


19.48(6) 
12 


20.3(2) 
5 


19.6(3) 
1.4 


18.5(6) 
0.18 
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Table 3: Fits of lattice data using 3-loop PT to i) g^-PT, ii) Leading Order 
Lattice Distorted PT, Hi) Next-to- Leading Order Lattice Distorted PT, and iv) 
9E-PT. 



8 





a 1 from 


A/309) 






M p 


u 


IP - IS 


gf^g^-Perturbation Theory 


Ai« [MeV] 
X 2 /dof 


111(1) 

0.483(6) 
4 


104(2) 

0.39(1) 

10 


69(6) 

0.04(8) 

1.4 


103(13) 
0.37(8) 
1.6 


120(30) 

0.47(14) 

0.1 


• 
• 


^^-Perturbation Theory 


A' v [MeV] 
d\ 

X 2 /dof 


• 
• 
• 


• 
• 
• 


110(10) 
0.05(5) 
1.4 


230(100) 

0.7(6) 

1.6 


• 
• 
• 


• 
• 


<7^^ -Perturbation Theory 


A' v [MeV] 
d\ 

X 2 /dof 


147(4) 

0.20(2) 

4 


122(3) 

0.01(1) 

14 


• 
• 
• 


120(20) 
0.01(8) 
1.6 


140(70) 

0.1(3) 

0.03 


0.15(1) 
4 


fl^-Perturbation Theory 


A E [MeV] 
df 

X 2 /dof 


• 
• 
• 


11.7(7) 
0.25(7) 
12 


• 
• 
• 


10(3) 

0.1(3) 

1.6 


• 
• 
• 


0.83(8) 
3 



Table 4: Fits of lattice data using eg. (177^) (with an unknown 3-loop parameter 
d 2 ) with g' defined as i) gjj s , H) gv' H \ an d Hi) gE2- A • signifies that no fit was 
found with < d 2 < 1. 
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6.5 



A CT [MeV] 5.5 - 



4.5 



~i i i i i i i r 



4 I I I I I I I I L 

5.4 5.6 5.8 6 6.2 6.4 6.6 6i 



Figure 1: Plot of A a = a a x fprido) against (3 = Q/gl- The 3-loop form of f PT 
was used. 



because there is not yet a great deal of data at a wide range of (3 values for 
non- Wilson actions. 



We also fit the discrete beta function A(3 defined implicitly by 

a L (p-Ap(J3)) = 2a L (J3f) 



(5) 



where is again defined from eq.(Q) with the same constraints on the coefficients 
as mentioned above. We use the A/3 data from the QCD-Taro collaboration [EH . 

In these fits we begin by including only the leading 0(a n ) term appropriate 
for each quantity. To be specific we rewrite eq.(Hf) as 



A, 



X 



x„ 



/?r(l) 



(6) 



where there is no implicit summation over n and v. Note that the 0(a n ) coeffi- 
cient has been normalised so that X n ^ is the fractional amount of the 0(gQa n ) 
correction at a standard value of go — 1, corresponding of course to /3 — 6/<7q = 6. 
Also note that we have approximated the polynomial coefficient, c^(^q), by its 
leading O(g^) term. For each physical quantity listed above this leading term 
corresponds to the following values of n and v: 



M P ,U,IP-IS 
a, A/3 
r 



v 
v - 



0,n = 1 

: 0,71 = 2 

2,71 = 2. 



(7) 



10 



6 



4 



2 



r r~ 

1 1 1 1 1 1 1 1 1 1 


- 




/ / 




/ 7 


Leading Order 


/ / 


Lattice Distorted-PT / 


/ 




/ 


~~ gM$bar" PT / / 


/ ■ 






/ / 




80 PT / 




/ ^ 




c "P^i i i 


i i 



5.5 6 6.5 7 

Figure 2: Plot of a^ 1 from the string tension together with fits from various 
methods described in the text. In all cases the 2-loop function fpr was used. 
The Monte Carlo data points are taken from table 0. (The next-to-leading order 
lattice- distorted PT curve is not shown, since it overlies the leading order lattice- 
distorted curve.) 



The values in eq.(|7]) arises because different quantities have different discretisation 
effects. For example, a, is measured using the gluonic part of the lattice action 
only, and hence is correct to 0(a 2 ). The scale r is correct to O^g^a 2 ) |§, hence 
the values v = 2, n = 2 for this quantity, etc. 

A simple least squares fit of each column of data in table [I] to eq.(H) provides 
the values for A^, the coefficient X n>u , and the \ 2 as listed in tables ^| & in the 
rows marked Leading-Order Lattice Distorted PT. Table |^ contains fits using the 
2-loop fpr-, and table ^ has fits using the 3-loop function.^ In the case of A/5, we 
use eq.(|]) with defined in eq.(|6]). 

Also shown in tables & are the fits to 2-loop and 3-loop go-PT.Q We see 
that leading order lattice-distorted PT fits the data very well compared to the 
g -PT case, with the \ 2 /dof down by up to two orders of magnitude. Figures 

4 Note that the fitting parameters for A/3 in table [2] correct a slight error in the corresponding 
values in ||. 

5 The functional form of this fit is simply a£ 1 (gg) = A^j /pt(<?o) (i.e. eq.(||) with X n _ v = 0). 
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o 

In 




5.5 6 6.5 7 





Figure 3: P/oi o/ a -1 /rom r together with fits from various methods described 
in the text. In all cases the 2-loop function fpr was used. The Monte Carlo data 
points are taken from table 



display the above fits in the case of the 2-loop fit. 

As a further check of the method, we include in the fit the next-to-leading 
term in a. However, due to the large statistical errors in some of the lattice data 
we perform this fit only for a, r and A(3 where the statistical errors are very 
small. The fitting function in these cases is: 



a 



^o 2 ) = >< 



/pt(0i 



r PT {i) n ' u+2 /pV 2 (i) 



The results of these fits are again displayed in tables & [5] for the 2-loop and 
3-loop cases. The figures [3] and [7| also display the fit in the 2-loop case for r and 
A{3. 

Obviously in the limit of infinite statistical precision, reasonable fits to the 
lattice-distorted PT formula would only be obtained if the 0(a n ) terms were 
included to all orders. The fact that it is necessary to go to next-to- leading order 
for the a, r and A{3 data to obtain a sensible x 2 /dof simply states that these 
data have sufficiently small statistical errors to warrant this order fit. From here 
on, we take the results of the next-to-leading order fits to the a, r and Aj3 data 
as our best fits for these quantities. 
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5.5 6 6.5 7 





Figure 4: P/oi o/ a -1 from M p together with fits from various methods described 
in the text. In all cases the 2-loop function fpr was used. The Monte Carlo data 
points are taken from table 



Some comments about the lattice-distorted PT procedure are necessary. 

• The most important point is that for the a, ro and A(3 data, the agreement 
between the data and lattice-distorted PT is truly remarkable considering 
the tiny statistical errors in the lattice data.0 

• Another important finding is that the values of for the various quantities 
in the 3-loop case are all consistent with Al = 7.7 MeV within la. The 
only exception is the string tension. Since the "experimental" value of 
the string tension requires certain model assumptions, it is not possible to 
draw any firm conclusion from this last observation. While the effects of 
quenching will mean that, at some level, the various values of A^ will not 
all coincide, it is encouraging that the A^ values agree to ~ 10% []. We 

6 In fact, a very close look at the second order fit to ro for the data from [jllj resolves a small 
discrepancy between fit and data for (3 = 5.9. It is possible that this is due to finite volume 
effects since L/a was 16 for (3 < 5.9 and 32 for j3 > 6.0 for this data fll5|| . 

7 Because the simulations were all quenched, it would be incorrect to perform a simultaneous 
fit to all the quantities {er, ro, M p , f v , IP — IS} using a single A^. 
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5.5 6 6.5 

P 



Figure 5: Plot of a -1 from f w together with fits from various methods described 
in the text. In all cases the 2-loop function fpx was used. The Monte Carlo data 
points are taken from table |7|. 
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5.5 6 6.5 

P 



Figure 6: Plot of a -1 from IP — IS splitting together with fits from various 
methods described in the text. In all cases the 2-loop function fpx was used. The 
Monte Carlo data points are taken from table 
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0.6 



0.5 



0.4 



0.3 




- — Next-to-Leading Order 
Lattice Distorted-PT 

Leading Order 

Lattice DiStorted-PT - 

gM0bar" PT 
g -PT 

J i I ! I I I I I I I I I L_ 



6.5 



7.5 



Figure 7: Plot of A/3(/3) together with fits from various methods described in the 
text. In all cases the 2-loop function fpx was used. The Monte Carlo data points 
are from [25]. 




0.05 0.1 0.15 0.2 0.25 0.3 

{a PT } 2 

Figure 8: Plot of a^ 1 /ap^ versus ap T for the string tension, where aJ, T = 
Kl/ fpT^gl) ■ The value of = 5.85 MeV from table [| was used and the 2- 
loop form of fpT was used to define apx- 
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take Al = 7.7± 10% MeV as an overall average where the error includes an 
estimate (obtained by a comparison of the 2 and 3-loop fits) of the effect of 
the unknown higher order terms b^, b\ etc. Using Ajjg = 28.81 x A^, [26], 
we have, 

A^T° = 220 ± 20 MeV, 

where the superscript Nf = refers to the quenched approximation. This 
compares very well with other lattice determinations of A^^° |57] and 
therefore is support for the validity of this approach. 



The typical values of X n>v in table ^| is 20-40%. In P8|, it was found 
that non-perturbative determinations of the renormalisation constant of 
the local vector current, Zy en , vary by 10-20% depending on the matrix 
element used. This spread in Zy en has been interpreted as 0(a) effects 



13 1 . Note that the simulation in p8[ was at /? = 6.4 where the value of a, 



and hence the size of the 0(a) term, is around half that at /3 = 6.0. Thus 
the values we obtain for X n>v in this study are directly comparable with the 
0(a) effects already uncovered in f28] , confirming the lattice-distorted PT 
picture. 

A compelling plot supporting lattice-distorted PT is shown in figure [S[ 
In this graph a~ x jap% is plotted against a 2 PT , where ajJp = Al/ /pt(<?o)- 
The 2-loop form of fpr is used, however, in this plot, the 3-loop would be 
indistinguishable. The value of A^ = 5.85 MeV from table |2] is used. If g$- 
PT were correct (i.e. if there were perturbative scaling) then the behaviour 
would be constant. Clear evidence for linear behaviour with non-zero slope 
is apparent, implying that the dominant corrections to go _ PT are of 0(ap T ). 

The coefficients for the second order terms, X n ^ u+ 2 are an order of magnitude 
smaller than the first order terms, X n ^ v . This follows our expectation that 
the expansion in fp-p in eq. @ forms a convergent series, and that the bulk 
of the cut-off effects at present values of (3 are due to the leading order 
term. 

One of the most exciting features of the lattice-distorted PT approach is 
that it can reproduce the behaviour of A/? (see fig.0)- 3-loop go-PT predicts 

A(3((3) = 6 + 6 l9 2 + b L 2 gt (9) 

with 5 = 126 log2, 8 X = 12&!log2 and 5% = Y2b\ log 2. The dotted 
curve in figure [7| with a near constant value of around 0.6, shows the be- 
haviour predicted from 2-loop go-PT. (The 3-loop term, 5 2 <?o brings this 
curve down by only around 2%.) The large discrepancy between lattice 



Monte Carlo data and 2-loop gcrPT for A/3 has long been noted [29|. The 



conventional reason for this "dip" is that it is a remnant of the line of first 
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order phase transitions in the fundamental-adjoint coupling plane which 
ends near the fundamental axis in the vicinity of this dip f30fl . However, 
the lattice-distorted PT can clearly solve this problem without resorting to 
other explanations. 

Note in p5] , a fit was performed to their Af3 data using two free parameters. 



The first was a 3- loop coefficient S2 (see eq.(gD), and the second involved 
defining a renormalised coupling, g u , with a shift c relative to the gjj§ 
coupling 0, 

6 6 
— — — c 

= ^%^ + 0.15 -c . 
% 

A good fit was obtained with a sensibly small 3-loop term, 62/ So = —0.013(1), 
but with a huge shift c = 2.61(4) which is unphysical since it leads to a 
divergent value of g 2 at (3 = 6/<7q ~ 5 PS[ . Thus this attempt at reconciling 
the A/3 data with PT proved unsuccessful. 

Comparing the values of X nv and X ntU+2 in the 2 and 3-loop cases (i.e. 
tables & 0) we see that they are compatible. In fact the main changes 
between the 2 and 3-loop cases are in the values of Ai which increases 
by around 15-20%. This variation is consistent with the value of d\ (see 
eq.(||)). In other words, the term d\g 2 can be well approximated by d\. 
This again implies that the higher order terms in g 2 are not important in 
the comparison between perturbation theory and current Monte Carlo data 
(apart from a normalising effect on A^). 

The values of X n>1/ and X n>u+2 in the case of a are compatible with those 
for Af3. This is a comforting feature, since both quantities are, in effect, 
obtained from the study of large Wilson loops [E5| . 



As far as the fit to M p , f n and the IP — IS splitting are concerned, the 
errors in the lattice data are large enough to allow many functional forms. 
Thus these data do not constrain the lattice-distorted PT fit (or fits from 
other schemes). 

The coefficients c' n of the 0(a n ) terms in eq.(£D are, strictly speaking, poly- 
nomial functions of g 2 . In the above fits, we have taken only the leading 
term in this polynomial. This is consistent with the lattice-distorted PT 
approach, since it assumes that the leading term of 0(a n ) appearing in 
eq.(|]) dominate the 0{g 2l ~ i ) terms. Consistency within this picture there- 
fore implies that the polynomials c^^q) can be replaced with their leading 
term in g^. 
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The errors in the fitting parameters X n>v are statistical only. They do not 
include any estimate of the systematic error due to the non-inclusion of the 
still higher order terms in a. A more refined fitting procedure would be 
required to estimate these errors. This will be left for future studies. 

The renormalisation constant used in the determination of aj^ is Z Ren {gfy = 
1 — 0.132(7q |3l] i.e. for consistency we have used the bare lattice coupling 



as the expansion parameter, rather than a renormalised coupling. However, 
as already mentioned, the statistical errors in the data for ajj are too large 
to enable useful conclusions to be drawn from these data alone. 



4 Fits Using a Renormalised Perturbation The- 
ory 

4.1 Introduction to the Renormalised Coupling fits 

In this section we fit the Monte Carlo data for a -1 to the following functional 
form: 

alHd) = T^WY (10) 

where g' is some "renormalised" coupling which is in turn a function of the bare 
coupling g . 

Note that the philosophy behind these fits is that the failure of asymptotic 
(perturbative) scaling is explained by higher order terms in perturbation theory. 
The coupling g' is defined with the aim of creating perturbative series with im- 
proved convergence properties. This is an orthogonal philosophy compared to 
the procedure outlined in Sect. |3| where finite lattice spacing errors are assumed 
to cause the disagreement between Monte Carlo data and perturbative scaling. 
(Note that the 0(a n ) terms in eqs.(|| & H) cannot be expressed as polynomials in 
<7o - i-e. they cannot be written in the form of the di_ l g'^~ A terms in eq.(f|).) 

The following subsections study fits using various definitions of g'. Both 2 
and 3-loop fits were performed, analogous to the lattice distorted case. At the end 
of this section we make some general comments about the success or otherwise 
of the renormalised coupling approach. 

4.2 Fits Using (/^-Perturbation Theory 



Following ref. J7J we define an MS— like coupling, where g' = gjjg, 

o ] , \ = ~2 < \ Tr Uplaq > MC +0.025. 
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The results of using this definition of g' in the fitting function Eq.flTOD (with 
the 2-loop definition of fpx) are displayed in Figs || - [7| and in Table || in the rows 
headed gj^. 

In order to perform 3-loop fits, an extra fitting parameter, d^f s , in fpr is 
required. This is because d\ from |l| is only appropriate for the go scheme. The 
results of these fits are displayed in Table [|. 



4.3 Fits Using gy-Perturbation Theory 

We now turn to some definitions of g' proposed by Lepage and Mackenzie M. 
Ref. 0] lists three alternative "renormalised" couplings gy (i = J, II, III) which 
are based on, or related to the strength of the static quark anti-quark potential. 
They are: 

4\ir/a) = - x (1 + 0.513a y ), 

< g 1 rUpi aq >MC 

(see § eq(29)); and 

-ln< ^TrU plaq > MC = y ^(3.41/a) x (1 - 1.19a y ), 
(see 0] eq(20)); and 

a{ / m) (46.08/a) = a L , 

(see M eq(30)) where a = g 2 /4n throughout and MC stands for Monte Carlo 
estimate. Lepage and Mackenzie argue that all three definitions agree up to 
0(a>y) (see discussion surrounding Fig. 7 of Ref. M). 

Before checking the fits of a -1 using these definitions of gv, we briefly discuss 
the effects of the momentum scale in the above definitions. Lepage and Mackenzie 
argue that there is a momentum scale q* which is most appropriate for each 
quantity studied which can be obtained from a mean field calculation. Thus, 
for example, the q* value for the critical quark mass for the Wilson action is 
q* = 2.58/a. In order to fit the Monte Carlo data for a -1 to Eq. ( |T0D (with 
g' = gy), we should, in theory, first obtain a value for q* appropriate for each of 
the quantities fl studied (Q = a,r , M p , f n , or the IP — IS splitting). However, 
as is shown below, fits using Eq. ([TDD are n °t dependent on the scale q* chosen, 
so long as the parameter Ay is trivially rescaled. 

We prove this statement as follows. Since gv is a function of q and j3, we write 
gv(q, ft) to make these dependencies explicit. In order to calculate the coupling, 
gv(q*,P), at some new scale, q*, we first must determine Ay. This can be done 
using 

Av = qfpT(g v (qif3))- (ii) 
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Once Ay has been determined, the value of gy can be determined at the new 
scale q* using the implicit definition: 

A v = q* f PT (g v (q\[3)). (12) 

Now we turn to the fit of the aj} obtained from a Monte Carlo determination 
of some physical quantity Q. The appropriate fitting function (see Eq . (|lOl) ) is 

a L (%) 



Using Eqs. ([□] & [TJ|) we immediately have 

! , _ q*/q Ay _ A' v (q) 
[90) fp T (9y(q,P)) fPT(9y(q,P)) ' 1 J 

where A' v (q) = —Ay. So the fit is independent of the choice of q* apart from 
a trivial rescaling of Ay, as claimed. Therefore we do not need to worry about 
the choice of q* for the quantities Q considered in this study. Note, however, 
that for quantities such as the plaquette which perturbatively are expressed as a 
polynomials in g 2 rather than proportional to J'pt, the above argument does not 
hold, and results do depend on the choice of q*. This last fact has been shown in 
Ref. Q. 

The results of the fits of the a -1 data to the functional form in Eq. ( |13|) (with 
the 2-loop definition of fpr) are shown in Table ||] in the rows headed gy' 1 ^ . 
The coupling gy H ^ was not considered since it would be the lead to the same fit 
as the go case apart from a trivial constant factor in the Ay value as discussed 
above. The q values for the fits I, II are q = ir/a, 3.41/a respectively. Since it 
is only the quality of the fits that is of interest no attempt was made to run the 
coupling gy' 11 ^ to the same momentum scale. However if this was attempted, 
inconsistent results would be obtained, since the A' v (q) values do not appear to 
fit the relationship A' y (g (/) ) = ^A' y (g (/7) ). 



As in the MS case, in order to perform 3-loop fits, an extra fitting parameter, 
d\, in f PT is required. The results of these fits are displayed in Table f|. 



4.4 Fits Using ^-Perturbation Theory 

Following refs. 0] and we define a coupling, g E , based on the plaquette, 



Cl 



g% I" < \T r Uplaq >MC 



where C\ — g is the coefficient of the first term in the perturbative expansion for 
1— < \TrU p i aq >. The results using this definition for g' in the fitting function 
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eg. (|T0|) are displayed in Table |2| in the rows headed g^-PT for the case of the 
2-loop fits. For the 3- loop fit, the coefficient appropriate for the gE scheme can 
be derived from a short calculation |l|], = 0.01161. The results for these fits 
are displayed in Table [| 

In ref ||, the gE scheme was extended to the next order in the perturbation 
series. That is, a coupling gE2, is defined by demanding that the Monte Carlo 
result for the plaquette is equal to the perturbative series truncated to second 
order, i.e. ^ 

1- < 7^TrU p i aq >MC= Cig E2 + C 2 g E 2, 

where c 2 = 2 x (0.0204277 - 1/288) ||. The results of this fit are displayed in 
Table |2] in the rows headed gE2 for the 2-loop fit. As in the MS case, in order 
to perform 3- loop fits, an extra fitting parameter, <if 2 , in fpp is required. The 
results of these fits are displayed in Table f|. 



4.5 Discussion of the Renormalised Coupling fits 

We begin with a discussion on the quality of the fits in the 2-loop case. As can 
be seen from Table |], the fits using gjjg, g v , Qy 1 , gE and gE2 -PT are not as 
successful as those from lattice distorted-PT. It is interesting to note that of the 
two gy fits, the g\J definition seems more successful. Surprisingly, the gE2 fit is 
not a great deal better than the qe fit, even though it should naively be correct 
to one more order of perturbation theory. This in itself suggests that pushing to 
higher orders in perturbation theory does not pay significant dividends. 

Concentrating on the fits to a, r and A/3 where the lattice data has very 
small statistical errors, we see from Table |^ that the gy and the gE2 -PT methods 
are the most successful of the renormalised coupling approaches. In fact the g v 
and gE2 -PT fits to r are comparable in quality to the leading-order lattice 
distorted PT. However, these same gy and gE2 fits are an order of magnitude 
worse for the a case compared to the leading-order lattice distorted PT, and are 
also significantly worse for the fits to Aj3. 

Of course the (leading order) lattice distorted PT fits include one extra fitting 
parameter (X n>v ) compared to the renormalised coupling fits, so it is tempting to 
argue that the former fits are more successful only for this reason. However, Sec.[5] 
details the results from a renormalised coupling style fit with two extra fitting 



parameters that proved entirely unsuccessful |25 . 

Turning to the 3-loop case, even more striking comments can be made. While 
it remains true that the lattice distorted PT fits remain stable and of high quality, 
the same cannot be said for the renormalised coupling approach. The gE scheme 
(whose 3-loop coefficient is known) cannot reproduce the quality of fits of the 
(even) the leading order lattice distorted PT (see Table 131). Furthermore, of the 
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four schemes listed in Table |], none can provide sensible fits for all the quantities 
studied. A "sensible fit" is defined, very liberally, as one where the fitted value 
for d,2 lies in the range < di < 1. Note also that in the cases where sensible 
fits are obtained, the values of d 2 depend strongly on the quantity being studied 
- a clear indication that there are non-universal (finite-lattice spacing) effects 
present. Furthermore, the values of the coefficient d 2 can often be ~ 20% or more. 
According to the renormalised coupling philosophy, d 2 should be significantly 
smaller than the corresponding coefficient in the g$ scheme. However, since d% ~ 
20%, this is clearly not the case. 

In conclusion none of the renormalised coupling methods approach the quality 
and consistency of the fits from lattice distorted PT method. 

5 Discussion & Conclusions 

In this paper we discuss the possible causes of the lack of perturbative scaling in 
dimensionful lattice quantities. We have argued that the cause of this disagree- 
ment is due to either (i) the neglect of higher order terms in perturbation theory 
(i.e. the renormalised coupling approach), and/or (ii) the presence of cut-off ef- 
fects due to the finiteness of the lattice spacing, a (i.e. the lattice distorted PT 
approach). Taking a collection of quenched lattice data from various collabora- 
tions for the string tension a, the hadronic scale r , M p , f n , the IP — IS splitting 
in charmonium and the discrete beta function Aj3, we have attempted fits assum- 
ing each of the hypotheses (i) and (ii) above. We have found that modelling the 
discrepancy using (i), through the use of a renormalised coupling, g', provides 
less satisfactory fits to the data compared to (ii). In fact, the hypothesis (ii) 
leads to a remarkable consistency both in the fitted A parameters and amongst 
the sizes of the 0(a n ) corrections. Furthermore, the hypothesis (ii) succinctly 
explains the "non-perturbative" behaviour of A/3. This, we feel, is one of its 
major achievements. Using the fits from lattice distorted PT (method (ii)) we 
estimate A^g^° = 220 ± 20 MeV, perfectly compatible with values obtained from 
other methods. We find excellent consistency for this value amongst the 
different physical quantities studied. 

It is clear however, that we have not proven that the lack of perturbative 
scaling is due to finite lattice spacing effects (i.e. hypothesis (ii)). In fact, pre- 
sumably, the true state of affairs is that the lack of perturbative scaling is due 
to a mixture of the two effects (i) and (ii). However, we have given strong ar- 
guments to support the claim that the discrepancy between the perturbative 
scaling formula and current Monte Carlo data is dominated by lattice artifacts 
(see Sec.0). 

A pictorial representation of the approach to perturbative scaling in lattice 
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Figure 9: Plot of aT 1 1 versus opt and g\. The string tension is used to set a. 
apT is defined as Apj /pt^o) where Al = 5.85 MeV is taken from Table [| and 
the 2-loop version of fpp is used. The data points are plotted using diamonds (the 
error bars are suppressed since they are smaller than the symbols). The vertical 
lines joining them to the apr — g% axis are guides for the eye only. The curves 
in the app — gl plane are app = fpT^gD/Ap with Ap defined as above. (No 
discernible difference in this plot would be observed if the 3-loop version of f PT 
were used instead.) 
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data is displayed in Figure [9|. Here the ratio a" 1 / ap^ is plotted on a 3-dimensional 
graph against both g% and apT- a' 1 is defined from the Monte Carlo values for the 
string tension (from Table [l]), and app is defined, Opr(<?o) = ^l/ fprido), where 
Al = 5.85 MeV is taken from Table @. (In the actual plot the 2-loop formula for 
fpT has been used. However, the 3-loop plot would be indistinguishable.) We use 
the string tension in this graph since the lattice data has very small statistical 
errors, however other quantities would lead to similar plots. The diamonds are 
the data points (error bars have been omitted because they are smaller than the 
symbols), and the vertical lines joining the data points to the app — g% plane 
are guides for the eye only. The two curves shown are defined by app^) = 
fpp(go)/A L , with A L defined as above. The lower curve is a guide for the eye, 
and the upper curve is in the a~ x jap\ = 1 plane. This upper curve shows 
the relationship between the Monte Carlo derived lattice spacing, a, the bare 
coupling, (?q, and the perturbative formula for the lattice spacing, apx in the 
idealised case of perfect (2-loop) perturbative scaling. It is important to note 
that the deviation between the data points and this curve in Figure |D| is precisely 
the discrepancy being studied in this paper. Q 

The ratio a~ 1 /api r is interesting since from eq.(f|) it is unity up to terms 
O(gl) and 0(a PT ), i.e. 

^ = (1 + E^-iSo'-T 1 x (1 + E c'M)a PT (9 2 ))- (14) 

a PT 1=4 n=l 

From Figure ^| the ratio a" 1 j(ip% varies from 0.6 to around 1.0 over the 
range of presently available lattice data. It can also be seen that the range 
in apT available from current lattice data spans a factor of 3-4, whereas the 
corresponding range in g$ is very small - from around 0.88 to 1.05 only (i.e. a 
spread of ^ 20%). Thus, the best of the current data is relatively a long way from 
the continuum limit a = g^ = in terms of variable g$, but apparently relatively 
close in terms of a. Therefore current lattice data is in the regime where effects 
of 0(a), 0(a 2 ) etc, will be seen clearly (if these terms are significant) since there 
is a broad range in a available. Conversely, higher order effects in g^ will be much 
more difficult to discern since the range available in this variable is so compact. 

We now aim to reconcile the rapid variation of a~ l /a~p\ from 0.6 to around 
1.0 with eq. ([l4]) . Considering the relative variations in g$ and apT, it is natural to 
assume that the 0(a n ) term (with n = 2 for the string tension) is more dominant 
than the O(g^) term. Note that the fact that significant 0(a) corrections have 
already been uncovered in lattice data adds support to these arguments (see 
Sec.|3|). 

8 Note that Figure [8] is a projection of Figure || onto the a^ 1 /op^ — apr plane. However the 
data in Figure @ is plotted versus a PT rather than apT- 
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To be more quantitative, we have performed a fit of the form 



i ?j\ 2, = (i + ^o 4 )" 1 x (i + 44tG?o 2 )), 

to the cr and r data in table |]. Obviously this fitting form allows both an 
0(a 2 ) distortion, and a 4-loop correction. The result is a worse fit than the next- 
to-leading order lattice distorted fit (which has the same number of free fitting 
parameters) in the case of the a data, and no sensible fit at all in the case of the 
ro data. 

Of course this argument again does not prove the case for lattice distorted 
PT; one could rely on the conspiracy of higher order terms in perturbation theory 
to recover the behaviour seen in Figure || without the need for any 0(a n ) terms 
at all. However, it seems much more natural, taking into account the arguments 
outlined above and in Sec.§, to rely on the 0(a n ) terms to explain the discrepancy 
instead. 

Using these assumptions, we can determine the value of the coupling when 
the finite lattice spacing effects first drop below, say, 1%. Taking a typical value 
ofX^ » 30% from Table | or [3], the 0(a n ) term in eq.(§, « X n f£ T (g%)/ f£ T (l), 
is less than 1% for g\ % 0.7, i.e. (3 ^ 9 in the case of n = 1 distorted quantities, 
and for g$ ^ 0.8, i.e. (3 ^ 8 for n = 2 quantities. 

Continuing the assumption that the discrepancy between lattice data and 
perturbative scaling is due to finite lattice spacing effects, we can also explain 
why it is that the renormalised coupling approach works as well as it does. This 
is presumably because the 0(a) effects of the quantity being used to define g' 
are embedded in the definition of g'. So long as the 0(a) effects in the quantity 
being used to define g' are similar in sign and magnitude to those in the quantity 
being studied (e.g. a, r , etc.) the PT expressed in terms of g' will fit the lattice 
data better than the bare go formula. Thus the g' definition mixes finite lattice 
spacing distortions with higher order perturbative effects. 

There is a final interesting opportunity assuming the correctness of the lattice 
distorted perturbative scaling philosophy |J . Suppose one is studying the contin- 
uum limit of a dimensionful quantity, fl, using lattice estimates determined at a 
set of finite a values. Typical examples of Q are the pseudoscalar decay constant 
f n or the mass M p . The physical value of Q can be written as 



n = lim \Z R -{gl) x n*{gl) x a^ 2 )] (15) 

3o->0 

where Z Ren is the renormalisation constant appropriate for Q (which is unity in 
the case of hadronic masses), and f2# is the (dimensionless) lattice value for Q. 
The limit go — * simply expresses the fact that the physical value is obtained 
by taking the continuum limit of the lattice estimates. Assuming that there is 
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no phase transition between values of gfjj used in present lattice calculations and 
g = 0, we can write eq.(|i"5|) as 



tt = \ lim Z 



^{gl)\ x lim p*{gl) x a^g*)} 

90^0 



For quantities where Z Ren {gfy can be expressed as polynomial in g^ (such as 
decay constants and hadronic masses) we have 



Q = Z Ren (g 2 = 0) x ^ 



where we have fitted a 1 according to eq.( 
along similar lines as follows, 



or eq. 



(16) 

and we have fitted fi* 



Ac 



x 



X„ 



PT\ 



Typically Z Ren {g\ = 0) is known immediately, without the need for any calcula- 
tion, (e.g. it is unity for decay constants and hadronic masses). This means that 
the continuum estimate of Q can be found using eq.fllBf) simply via two lattice 
distorted PT style fits. 
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